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Gigantic enhancement of spin Seebeck effect by phonon drag 
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We investigate both theoretically and experimentally a gigantic enhancement of the spin Seebeck effect in a 
prototypical magnet LaY2Fe50i2 at low temperatures. Our theoretical analysis sheds light on the important 
role of phonons; the spin Seebeck effect is enormously enhanced by nonequilibrium phonons that drag the 
low-lying spin excitations. We further argue that this scenario gives a clue to understand the observation 
of the spin Seebeck effect that is unaccompanied by a global spin current, and predict that the substrate 
condition affects the observed signal. 
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When a temperature gradient is applied to a ferromag- 
net, a force is induced acting on electrons' spin to drive 
spin currents. This phenomenon termed the spin Seebeck 
effect (SSE)i"2 has recently drawn tremendous attention 
as a new source of spin currents needed for future spin- 
based electronics^^ SSE is now established as an uni- 
versal aspect of ferromagnetic materials as it has been 
observed in a variety of materials ranging from a metal- 
lic ferromagnet 1 NisiFeig and a semiconducting ferro- 
magnet 2 GaMnAs to an insulating magnet 3 LaY2FesOi2. 
Besides its impact on the technological application, SSE 
offers a number of new topics on the interplay of heat 
and spin currents and it triggered the emergence of 
the new field named "spin caloritronics" 8 in the rapidly- 
growing spintronics community. 

A mystery concerning SSE was how conduction elec- 
trons can sustain the spin voltage over so long range of 
several millimeters in spite of the short conduction elec- 
trons' spin- flip diffusion length A s f, which is typically of 
several tens nanometers. This problem has recently been 
resolved by a series of experiments on spin currents using 
magnetic insulators. A recent experiment on the elec- 
tric signal transmission through a magnetic insulator 9 
highlights the role of the low-lying magnetic excitation 
of localized spins, i.e., magnons, by demonstrating that 
magnons transmit the spin current over a long distance 
of several millimeters. A subsequent experiment on SSE 
for a magnetic insulator 3 LaY2FesOi2 confirmed that the 
magnon-based scenario can explain the SSE experiment 
at room temperature, since the length scale associated 
with magnons 3> A s f. However, a new issue on SSE was 
brought by a very recent experiment on a ferromagnetic 
semiconductor 2 GaMnAs, where it was demonstrated, by 
cutting the magnetic coupling in GaMnAs while keep- 
ing the thermal contact, that SSE can be observed even 



in the absence of global spin current flowing through 
GaMnAs. Obviously, the scenario of magnon-mediated 
SSE 10 ! 11 fails to explain the experiment, showing that 
the full understanding of SSE has not yet been reached. 

For a deep understanding of the physics behind SSE, 
we here explore the low-temperature behavior of SSE in 
an insulating magnet LaY2Fe5 0i2- Figure Q] shows a 
schematic illustration of our device structure^ 2 - An in- 
plane external magnetic field H and a uniform tempera- 
ture gradient VT were applied along the z direction [see 
FIG. [1] (a)]. The VT generates a spin voltage across the 
LaY 2 FesOi2/Pt interface, and injects (ejects) a spin cur- 
rent 7 S into (from) the Pt wire. In the Pt wire, a part of 
the injected/ejected I s is converted into a charge voltage 
through the so-called inverse spin-Hall effect (ISHE): 13 



Vishe = e„{\e\I s )(p/w), 



(1) 



where |e|, Qh, P and w are the absolute value of electron 
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FIG. 1. (Color online) Gigantic enhancement of SSE in 
LaY2FesOi2 at low temperatures, (a) Schematic illustra- 
tion of the LaY2FesOi2/Pt sample and the temperature pro- 
file along the z direction. Here H denotes an external 
magnetic field (with magnitude H). The sample comprises 
a LaY2FesOi2 film with 8 x 4 mm 2 rectangular shape and 
two separated Pt wires with the width w attached to the 
LaY2FesOi2 surface at the interval of 5.6 mm. (b) T depen- 
dence of V/AT at H = 100 Oe. 
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FIG. 2. (Color online) Diagrammatic representation of the 
thermal spin injection process, (a) Magnon-mediated SSE. 
Here the system is composed of ferromagnet (F, in the ex- 
periment LaY2FesOi2) and nonmagnetic metals (N, in the 
experiment Pt), which are divided into three temperature do- 
mains of F\/N\, F2/N2, and F3/N3 with their local temper- 
atures of Ti, T, and T3. The thin solid lines with arrows 
(bold lines without arrows) are electron (magnon) propaga- 
tors. Here, J st j (J cx ) is the strength of the s-d coupling at 
the F/N interface (the exchange coupling in F). (b) Phonon- 
dragged SSE where the dashed lines are phonon propagators. 
The process P4 injects the spin current with the same mag- 
nitude as (but opposite sign to) the process P2 due to the 
relation T — T2 — — (T3 — Ta), while no spin current is in- 
jected into N2 because of the cancellation b etween the two 
relevant processes P3 and P3. Here, Q.o = y/ K p h/Mi OIl with 
the ion mass M[ on and the elastic constant K p h in F. (c) Cal- 
culated spatial dependence of the spin current injected into 
Ni (i = 1, 2, 3). 



charge, spin-Hall angle, resistivity and width of the Pt 
wire, respectively. Therefore, the VT-driven spin injec- 
tion, or SSE, is electrically detectable. In FIG. Q] (b), 
we show the temperature (T) dependence of Vishe/AT 
at H = 100 Oe, measured when the Pt wires are at- 
tached to the lower- and higher-temperature ends of the 
LaY 2 Fe 5 0i2 layer, respectively. The sign of V/AT is 
reversed between these Pt wires, a situation consistent 
with SSE- induced ISHE. Notable is that, at T = 50 K, 
the magnitude of V/AT is dramatically enhanced. 

A simple scenario of the magnon-mediated SS E 10 i n 
[FIG. (a)] is unable to explain the observed low-T 
enhancement. If such a scenario could explain the ex- 
periment, the low-T enhancement of V[she would come 
from cither the enhancement of the spin-Hall angle 0# or 
that of the magnon lifetime; otherwise I s due to magnon- 
mediated SSE is a monotonic increasing function of T as 
discussed below. From the T-dependence of the spin-Hall 
conductivity^ we conclude that there is no enhancement 
of Qh at low T. While the possibility of the enhance- 
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FIG. 3. (Color online) Comparison of experimental and the- 
oretical SSE signal. Solid circles: experimental spin See- 
beck data for LaY2FesOi2 (Dotted line is a guide to the 
eye). The solid curve: calculated T-dependence of Vishb 
due to the sum of the phonon-dragged SSE and the magnon- 
mediated SSE. The dashed curve: calculated T-dependence 
of Vishe due only to the magnon-mediated SSE. We have as- 
sumed T-independent Qh and a, and used Td = 565 K— and 
Ta/ = 560 K. The data are normalized by its value at 50 K 
except that the result for magnon-mediated SSE is plotted to 
reproduce the room-temperature signal. Inset: experimental 
thermal conductivity k for Y3FesOi2 taken from Ref.(2^ (solid 
circles) and the result of the fit (solid curve) using Eq. Q . 



ment of magnon lifetime is not conclusively excluded, 
judging from the ferromagnetic resonance linewidth in 
YsFesOi^ as a measure of the inverse magnon lifetime, 
it does not seem to be the case. Therefore, we need a new 
mechanism to account for the observed low-T enhance- 
ment of SSE. 

Here, the gigantic enhancement of SSE for LaY2Fe50i2 
below room temperature is analyzed in the light of 
phonon-drag mechanismji£~— Back in 1946 in the con- 
text of thermoelectrics, Gurevich pointed outiS that the 
thermopower can be generated by a stream of phonons 
driven by the temperature gradient, which then drag 
electrons and cause their convection. This idea, known 
nowadays as phonon-drag mechanism, has been estab- 
lished— iiS as a principal mechanism causing the low-T 
enhancement 19 of the thermopower. Because SSE is a 
spin counterpart of the Seebeck effect, it is natural to 
expect that a similar physics underlies SSE. It is this 
approach that we adopt in the present work. 

Our theoretical analysis starts from considering the 
model shown in FIG. (a). The key point in our model 
is that the temperature gradient VT is applied over the 
ferromagnet, but there is locally no temperature differ- 
ence between the ferromagnet (F) and the attached non- 
magnetic metals (N) . We assume that each temperature 
domain is initially in local thermal equilibrium, then we 
switch on the interactions among the domains and calcu- 
late the nonequilibrium dynamics of spin density in N. 

The central quantity that characterizes SSE is the spin 
current I s injected into N (in experiment Pt), since it 
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FIG. 4. (Color online) Schematic illustration of SSE unac- 
companied by a global spin current. The phonon-drag pro- 
cess which explains the experiment^ is shown. The meaning 
of each line (propagator) is the same as in FIG. [2] 



is proportional to the experimentally-detectable electric 
voltage via ISHE [Eq. Q}]. Following Ref. [TJ the spin 
current /™ as (iVi) injected into N\ due to the magnon- 
mediated SSE^^ is calculated as 

(2) 

where a is the Gilbert damping constant, Tm is the 
characteristic temperature corresponding to the magnon 
high-energy cutoff, and P is a nearly-temperature- 
independent coefficient. 20 Equation ^ means that the 
magnon-mediated SSE cannot explain the low-T en- 
hancement of the signal (the dashed curve in FIG. 

Now we proceed to a detailed analysis of the In- 
dependence of SSE in terms of the phonon-drag mech- 
anism. The Feynman diagram for the phonon-drag pro- 
cess in the present situation is shown in FIG.[3](b), where 
the phonons feel the temperature difference between Fi 
and F 2 , and drag magnons through the magnon-phonon 
interaction. Since the nonequilibrium phonons affect the 
magnon dynamics, this process injects spin current into 
Ni . The important point is that the spin current /^ rag in- 
jected in this process becomes proportional to the phonon 
lifetime r p h as^ 

lf^(N 1 )/AT = P , T ph B 1 B 2 , (3) 
where Si = {T/ J?/ ^ D,T du . j". .„, and B 2 = 

1 4ir J JO smh' ! (ii/2) £ 

^(^) 3 f /T ^ with the Debye 
temperature To, and P' is a nearly-temperature- 
independent coefficient^ Since r p h in a high-purity spec- 
imen is known to increase steeply at low T because of the 
rapid suppression of umklapp scattering^ it leads to the 
drastic enhancement of the phonon-dragged SSE. In our 
analysis, the T-dependence of r p h is extracted from the 
thermal conductivity data for YsFesOi? 2 ^ (see the inset 
of FIG. O using 2 !^ 

k(T) = (l/3)v 2 ph C ph (T)r ph (T), (4) 



where i; p h is the phonon velocity, and C p h(T) — 
9N D k B (T/T D ) 3 £ D/T dw isi J 2 \ w/2) is the phonon spe- 
cific heat with the number of phonon modes Njj . After 
getting the information on r p h(T), we calculate the In- 
dependence of T^ishe resulting from the phonon-dragged 
SSE. The result, plotted in FIG. [3] (the solid curve), 
shows an excellent description of the low-T enhancement 
of SSE. 25 Our analysis demonstrates that the phonon- 
drag mechanism is of crucial importance to understand 
SSE below the room temperature. 

Finally, we show in FIG.Hlour interpretation on the ob- 
servation of SSE that is unaccompanied by a global spin 
current, 2 where the heat is carried by phonons through 
the nonmagnetic substrate while the spin is injected lo- 
cally at the F/N interface. This interpretation is rein- 
forced when we recall that the magnitude of the spin See- 
beck signal is enhanced with decreasing T even well be- 
low the Curie temperature, whose tendency is consistent 
with the phonon-drag mechanism as is seen in FIG. [3] 
Furthermore, the fact that the experiment was done be- 
low the room temperature supports the the phonon-drag- 
based scenario, since the phonon-drag process becomes 
more effective at low T as emphasized in the previous 
paragraph. All these considerations strongly support 
that the SSE experiment for GaMnAs can be interpreted 
in terms of the phonon-drag mechanism, and results in 
a prediction that the substrate condition affects the ob- 
served signal. Our demonstration opens a new route to 
control spin currents by means of phonons and stimulates 
further progresses in spin caloritronics. 
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1. Experimental details 



2. Derivation of Eq. (3) 

Following Ref. [SI], the spin current I s (Ni,t) injected 
into the nonmagnetic metal Ni (i = 1, 2, 3) is calculated 
as 



I s (Ni,t) = - 



4j: 



k — q 



ft y/2N F N N h' 



ReC< q (t,t), 



(SI) 



where Np (Nn) is the number of lattice sites in F (N), 
Sq is the size of the localized spins in F, and J^ q is 
the Fourier transform of the s-d interaction at the F/N 
interface. Here, C^ q (t,t r ) — — i(a q (t')s^ (*)) measures 
the correlation between the magnon operator a q and the 
spin-density operator s^ = (s% — is v k )/2. Note that the 
time dependence of I s {Ni,t) vanishes in the steady state 
and it is hereafter discarded. Introducing the frequency 
representation C< q (t-t>) = &C< fc Me-M*-* ) 

and adopting the representation [S2] C = I ' CA J as 

C R + C A ], we 



well as using the relation C 
obtain 



< -\\c 



K 



q.h 



k — q 



So 



V>2N F N N h 



2^ 



(S2) 



for the spin current I s (Ni) in the steady state. 

When we introduce a renormalized magnon propagator 
5X q (u)), the interface correlation C appearing in Eq. (fS2"j) 
is generally expressed as 



q k-q 



y/N N N F fl 



Xk{u)6X q (u), 



(S3) 



where x k (u) = 



is the spin-density propa- 



gator satisfying the local equilibrium condition: 



x£M = 2iIm X £Hcoth(^) 



2k B T' 



(S4) 

Here the retarded component of Xfe( w ) i s given by 
Xfc( w ) = Xiv/(1 + X^k 2 - iujT si ) [S3] where xn, Kf, and 
r s f are the paramagnetic susceptibility, the spin diffusion 
length, and spin relaxation time. 

We now consider the phonon-dragged SSE [the pro- 
cess Pi shown in FIG. 2 (b)]. The renormalized magnon 
propagator 5X q (uj) in the present case is given by 



The single-crystal LaY2FesOi2 (111) film with the 
thickness of 3.9 /zm was grown on a GdaGasO^ (HI) 
substrate by liquid phase epitaxy. The 15-nm-thick Pt 
wires were then sputtered in an Ar atmosphere. The 
length and width of each Pt wire are 4 mm and 0.1 mm, 
respectively. The temperatures of the lower- and higher- 
temperature ends of the sample were respectively stabi- 
lized to T and T + AT, where T was controlled in the 
range of 300-50 K by means of a closed-cycle helium re- 
frigerator. 



6X q (u) = X q (w)t q {u)X q (u), 



(S5) 



o x A (J) J 1S magnon prop- 

agator satisfying the equilibrium condition: 



where X q (oj 



2ilmX^(uj)coth(^ T ). 

(S6) 

Here, the retarded component is given by X q (u) — 
(lu — u) q + iauj)^ 1 , where ui q = ^Hq + tu q is the magnon 
frequency for uniform mode 7-f/o a- n d exchange mode 
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D CK q 2 /h. In Eq. ([55]). the selfenergy S due to 



phonons is given by 



Eq(w) 



— y 



2N F ^ \ h 

K 



5D R (v)X q _ 



+ ^ (v)fi*,_ (w_ ) + 5Z? A ' (w_ ) } , (S7) 



where ^K,q = 



guj q 



Hut, 



-5— is the magnon-phonon in- 



teraction vertex with vk, v p h and M; on being the phonon 
frequency, phonon velocity and the ion mass, f is the 
Pauli matrix in the Keldysh space, and we have intro- 
duced the shorthand notations w_ = cj — v, q = q — K, 
and J = In Eq. (jS7|) . the full phonon propaga- 

tor SDk [the whole of the phonon lines for P2 in FIG. 2 
(b)] is written as [S4] 



8D K (u)=5D^ 9 (u) + 6D^(u). 
Here, SD l - e i{uj) = 



(S8) 



'SD'-^ R ( V ),6D'-">- K {u)\ ■ , 

0, 5D l ->«< A (v) J lh LUL 

local-equilibrium propagator satisfying the local- 
equilibrium conditions 5D L £ q ' A {v) = [5D l £ q ' H {v)]* and 
5D l ^ K {u) = [6D l K ^ R (u) - «5D^(^coth(^) 
with its retarded component given by 



the bare phonon propagator satisfying the equilibrium 
condition: 

= P&Mh = 2iImD« (1/) coth(^), 

with its retarded component and the phonon lifetime 
given by D^v) = (^-^ + i/r ph ) _1 -(^+^ + i/Tp h ) _1 
and r ph . 

Now substituting these expressions into equation (|S5j) 
and use Eq. (jS2|) . the spin current injected into Ni by 
the phonon-drag process is calculated as 



R 



N N N F 



y / dvKV^iTK^) 2 A k . q (v) 

1„ „ Ts" <J 



x[ coth (fe)-coth(^)], (S12) 

where i? = J^5o)^AWas/A)T ph /(4^ 3 fi 4 ^) with 
vd = Vph/ a s , and 



A M (i/) = / lm X H ^)lmX R _{^)\X R {^ 

x [coth(^) - coth(^)]. (S13) 



By setting T = T 2 , AT = T± — T 2 , Qq = v D and after 
some algebra, we obtain Eq. (3) in the main text. 
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N F (A/a s ) 



[ D kH] 2d k'(v), (S9) References 



if' 



while 8D n ~ eq — ( jj' 5£> ) is the nonequilibrium prop- 
agator with its Keldysh component given by 
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n-eq,K 
K 
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N F (A/a s ) 
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